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Abstract. As suggested by some extensions of the Standard Model of particle 
physics, dark matter may be a super- weakly interacting lightest stable particle, while 
the next-to-lightest particle (NLP) is charged and meta-stable. One could test such a 
possibility with neutrino telescopes, by detecting the charged NLPs produced in high- 
energy neutrino collisions with Earth matter. We study the production of charged 
NLPs by both atmospheric and astrophysical neutrinos; only the latter, which is 
largely uncertain and has not been detected yet, was the focus of previous studies. We 
compute the resulting fluxes of the charged NLPs, compare those of different origins, 
and analyze the dependence on the underlying particle physics setup. We point out that 
even if the astrophysical neutrino flux is very small, atmospheric neutrinos, especially 
those from the prompt decay of charmed mesons, may provide a detectable flux of 
NLP pairs at neutrino telescopes such as IceCube. We also comment on the flux of 
charged NLPs expected from proton-nucleon collisions, and show that, for theoretically 
motivated and phenomenologically viable models, it is typically sub-dominant and 
below detectable rates. 
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1. Introduction 

The fundamental nature of dark matter poses a profound challenge to contemporary 
theoretical particle physics. Observations constrain the neutrino — the only electrically 
and color neutral non-baryonic elementary particle within the Standard Model of particle 
physics — to have a negligible contribution to the overall dark matter budget [1,2]. Dark 
matter is regarded as one of the most compelling hints towards new physics beyond 
the Standard Model. The question of its elementary essence has triggered enormous 
theoretical and phenomenological efforts [3]. 

The existence of suitable dark matter particle candidates in several theoretically 
cogent extensions of the Standard Model, like low-energy supersymmetry [4] or extra- 
dimensional scenarios [5, 6], focused a strong interest on a class of dark matter candidates 
known as weakly interacting massive particles (WIMPs). Similarly to other Standard 
Model particles, WIMPs would fall out of thermal equilibrium and freeze out in the 
early Universe, leaving a relic abundance compatible with the inferred amount of dark 
matter [7]. These WIMPs can be directly detected by experiments looking for the 
minuscule energy deposition caused by dark matter particles scattering nuclei [8]. The 
pair annihilation of WIMPs into energetic gamma rays, neutrinos, and antimatter, is a 
second, yet indirect, handle on the presence and potential imprint of galactic particle 
dark matter [9]. 

The connection of the aforementioned scenarios to the electroweak scale, soon 
to be probed at the Large Hadron Collider (LHC), motivated the exploration of 
complementarity between collider searches for new physics and the question of the 
elementary nature of dark matter [10-14]. In most cases, if dark matter is a WIMP, 
the anticipated experimental signature at LHC would be the production of strongly 
interacting massive particles which promptly decay to the lightest and stable WIMP, 
plus a number of energetic jets and leptons. The neutral dark matter particle 
would escape the detector unobserved, leading to large missing transverse energy as 
well [15]. Conclusive identification of escaping neutral particles at LHC with dark 
matter permeating our Galaxy and other cosmic structures would, however, require 
some evidence from the direct and/or indirect WIMP searches listed before [9]. 

WIMPs are indeed attractive dark matter candidates, but are not the only 
theoretically envisioned possibility. The dark matter particle could exhibit even feebler 
interactions with ordinary Standard Model particles than a WIMP, making direct and 
indirect searches completely hopeless. For instance, the supersymmetric gravitino [16] 
or the Kaluza-Klein graviton of universal extra dimensions [6,17] are perfectly plausible 
'super-weakly interacting' [18] dark matter candidates (super- WIMPs). If Nature chose 
an option like this, collider signatures of new physics, if any, would strongly depend upon 
the nature of the next-to- lightest supersymmetric particle (NLSP). Since a super- WIMP 
is also very weakly coupled to the other new-physics heavier states, the NLSP would 
likely be quasi-stable. If the NLSP is neutral, the qualitative experimental landscape 
would look like that of a standard WIMP scenario. However, we would lack the needed 
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proof of a connection between the weakly interacting long-lived particles produced at 
colliders and galactic dark matter. 

If the NLSP is instead charged (constituting a charged massive particle, or 
CHAMP), the LHC would potentially observe the extremely distinct signature of a 
'heavy muon': charged tracks and penetration of the outer muon sub-detector, possibly 
at very low relativistic beta. CHAMPs are constrained by direct collider searches at 
LEP2 [19], as well as at the Tevatron [20]. \ 

If a CHAMP were stable on collider scales, it could still decay on cosmological 
scales, and thus impact precision astrophysical measurements [22], including the 
chemical potential associated with the cosmic microwave background black-body 
spectrum [23], the extragalactic gamma-ray background [24], the reionization history 
of the universe [25], the formation of small scale structures [26,27] and the synthesis 
of light elements in the early Universe [28, 29] (see, for implications of neutral particle, 
Ref. [30]). Anomalies in the above mentioned quantities, however, could hardly be 
considered smoking gun evidence that the collider CHAMPs were indeed related to 
lighter, super-weakly interacting dark matter. 

If CHAMPs were stable on collider timescales, but featured a short lifetime on 
cosmological timescales, say on the order of a year or less, CHAMPs produced in colliders 
might be trapped in large water tanks surrounding the detectors [31]. The tanks would 
then be periodically drained to underground reservoirs where CHAMP decays might 
be observed in low-background conditions [31]. While certainly not straightforward 
experimentally, such a technique might provide important information first on the actual 
met a- stability of the charged species, and, secondly on the nature of the super-weakly- 
interacting particle the CHAMP would decay into. However, even if CHAMP decays 
were actually observed, this would still not suffice as conclusive evidence that the elusive 
particle CHAMPs decay into is indeed the dark matter constituent. 

To our knowledge, beyond high-energy collider experiments, the only direct 
experimental handle on a super-weakly interacting dark matter particle featuring a 
heavier, meta-stable charged partner is CHAMP pair production via neutrino-nucleon 
collisions, followed by direct observation at neutrino telescopes. This idea, originally 
proposed in Ref. [32] , relies on the fact that the energy losses of CHAMPs in Earth are 
significantly smaller than those of muons, therefore CHAMP pairs (unlike muon pairs) 
can reach the detector even if they were produced far away. This makes the relevant 
target volume for neutrino-nucleon interactions much larger. The CHAMP pairs can be 
efficiently separated from muon pairs, due to large track separations in the detector. The 
original proposal was subsequently followed up by related studies [33-38], which focused 
on the specific case of a gravitino lightest supersymmetric particle (LSP), and a stau 
NLSP playing the role of the CHAMP. Among other aspects, these studies investigated 
in detail stau energy losses in Earth and in the detector, computed expected event rates 
for a few sample models and the relevant background, and addressed the possibility of 

| These and future searches at LHC are not trivial, but advanced work is being done to ensure that 
such potential signals would not be missed [21]. Detection for masses < 1 TeV is essentially guaranteed. 
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discriminating single-stau from single-muon events. 

Given the steady progress in the deployment of next-generation km 3 -size neutrino 
telescopes — particularly IceCube at the South Pole, already under construction and 
taking data — we consider it timely to address in detail a few points relevant for 
improving our understanding of the prospects for implementing the above outlined 
technique. In particular, as background rejection is not a substantial issue, the crucial 
point appears to be the evaluation of the CHAMP pair event rate at IceCube. To this 
end, we focus on the following four aspects: 

(i) Incoming neutrino flux: So far, all long-lived CHAMP analyses for neutrino 
telescopes have considered a relatively optimistic flux of astrophysical neutrinos — as 
large as the Waxman-Bahcall (WB) bound [39] — as the primary source for CHAMP 
pair production. However, the WB bound is only a theoretical upper limit on 
the flux of astrophysical neutrinos expected from optically thin sources; therefore, 
the absolute normalization as well as the spectral shape of the true astrophysical 
neutrino flux remain largely unknown (of course the biggest reason for this is that 
these neutrinos are as yet undetected!). On the other hand, atmospheric neutrinos 
have been detected and their flux is rather accurately known (below ~100 TeV [40]). 
At larger energies (> 100 TeV), although there is no detection so far, one can 
rather reliably extrapolate the flux of conventional atmospheric neutrinos from 
lower energies. In addition, one also expects a significant flux of so-called prompt- 
decay atmospheric neutrinos, which originate from the decay of short-lived charmed 
mesons and feature a harder spectral index than the conventional component. While 
we again have only an upper limit on this prompt decay component, we know 
from particle physics that it necessarily exists, at some level, in the high-energy 
regime, and will be measured accurately by IceCube. In any event, atmospheric 
conventional and prompt-decay neutrinos evidently contribute as well to CHAMP 
pair production in neutrino-nucleon collisions. In this paper, we study the role of 
these standard, guaranteed neutrino sources, and compare it to the contribution 
from astrophysical neutrino flux models including the WB bound. 

(ii) Underlying particle physics model: The event rate depends not only on the incoming 
neutrino flux, but also on the nature of the assumed particle physics model. Here, 
we consider generic supersymmetric models featuring a gravitino LSP and stau 
NLSP, and study how the stau pair production cross section and event rates 
at neutrino telescopes depend on the given mass spectrum. Neutrino-nucleon 
interactions produce slepton-squark pair final states, the amplitudes mediated by 
supersymmetric fermion exchange (neutral and charged gauginos and higgsinos). 

(iii) Other CHAMP sources: Unlike the high-energy neutrino flux, the flux of very 
energetic protons is accurately measured and well-known up to extremely high 
energies; proton-nucleon interactions feature a supersymmetric pair production 
cross section significantly larger than that of neutrino-nucleon processes. It 
therefore seems reasonable to quantitatively assess the flux of stau pairs produced 
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in the interaction of incident primary protons with nuclei in the atmosphere. The 
trade-off is the enormous cross section for proton-nucleus scattering into Standard 
Model particles, which depletes the incoming proton flux, highly suppressing 
any stau pair production rate. Since quark-ant i- quark processes can directly 
produce stau pairs, however, the kinematic threshold for stau pair production, 
as a function of the incoming primary particle energy, is lower for proton- 
nucleus than for neutrino-nucleus processes. Yet we find that, for reasonable and 
phenomenologically acceptable particle models, the expected stau pair event rates 
from proton-nucleus collisions are experimentally negligible at all energies, leaving 
neutrinos as the only relevant primary source particles for CHAMP pair production. 

(iv) Simplified analytic approach: We present the computation of the flux of staus from 
neutrino-nucleon interactions from first principles, and argue that, to an acceptable 
degree of accuracy, the total number of expected staus can be computed as one 
simple integral of three factors. Specifically, we show that the quantities of physical 
relevance are (1) the incident flux of primary neutrinos, (2) the ratio of the cross 
sections of neutrino plus nucleon into supersymmetric particle pairs to the total 
neutrino-nucleon cross section, and (3) a geometric efficiency factor. 

Hereafter, we specifically use supersymmetric staus as charged meta-stable NLPs, 
but note that the following arguments are applicable to any other possible candidates 
of long-lived CHAMPs. The outline of the remainder of the paper is as follows. We 
discuss the various components of the high energy neutrino flux and their connected 
uncertainties in Section 2. We introduce the new physics scenarios and compute the 
stau pair production cross sections in Section 3. Section 4 outlines the computation 
of the stau event rate at IceCube, including the above-mentioned simplified analytic 
treatment. The stau rate dependence on the particle physics framework is addressed in 
Section 5, and we draw conclusions in Section 6. 

2. High-energy neutrino flux 

In this Section, we summarize the high-energy neutrino fluxes we consider in the present 
study. As discussed in Section 1, past works considered only a flux of neutrinos close 
to, or saturating, the WB upper limit. However, other neutrino sources potentially 
contribute as well: these include conventional atmospheric neutrinos, prompt-decay 
atmospheric neutrinos and possibly astrophysical neutrinos other than those considered 
in the WB setup. Since the stau production rate does not depend on neutrino flavor, § 
all we care is the total neutrino plus anti-neutrino flux, i.e., the flux of v = v e + + z/ T , 
where each z/j here indicates neutrino plus anti-neutrino of flavor i. From this point on, 
we mean this combined quantity when we use the term 'flux', unless otherwise stated. 

§ Flavor is conserved in the underlying supersymmetric particle pair production event, but all 
supersymmetry particles then decay promptly to a stau plus Standard Model particles. Thus, an 
electron neutrino eventually produces a stau pair just as a tau neutrino does. 
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In Figure 1, we summarize and collect the various neutrino sources we discuss below, 
and the ranges of normalizations we consider. 

2.1. High-energy astrophysical neutrinos from extra- galactic sources 

Very powerful astrophysical objects such as active galactic nuclei (AGNs) and gamma- 
ray bursts (GRBs) are candidate high-energy neutrino sources. This is because strong 
gamma-ray emission detected from these objects can be attributed to the particle 
acceleration and successive interaction with the surrounding medium, magnetic and 
photon fields, which might also be a source of neutrinos via charged meson production. 

If these neutrino sources are optically thin, then the upper bound on neutrino flux 
is obtained from the well-measured cosmic-ray flux, because each proton that arrives 
at Earth should produce no more than a few neutrinos at the source. Based on this 
argument and assuming that the cosmic-ray spectrum above 10 9 GeV is of extra-galactic 
origin, Waxman and Bahcall [39] derived the upper bound for the flux before neutrino 
oscillation to be El d<&^ B jdE Vfi = (1-4) x 1CT 8 GeV cm~ 2 s _1 sr _1 , where the range 
reflects cosmological evolution of source density. As we expect that flavor ratio at 
production (i.e., before oscillation) is v e :i' ll :i/ T = 1:2:0 due to meson decays, the WB 
bound summed over flavors after neutrino oscillation is E^d&^ B /dE u = (1.5-6) x 10~ 8 
GeV cm" 2 s" 1 sr" 1 . Here we adopt E 2 v d^ B /dE v = 5 x 1(T 8 GeV cm" 2 s" 1 sr" 1 for 
our reference value and show this bound in Figure 1. 

We stress that although all previous studies adopted the WB upper limit on 
astrophysical neutrino flux as their incident source, that flux is not a model prediction, 
but an upper bound. One can never regard the output of a WB maximal neutrino 
flux as a solid prediction for the stau event rate at IceCube; the resulting stau flux is 
merely an upper limit. In addition, a few comments are in order on the WB bound's 
robustness, especially at our energies of interest. First, the WB bound is valid for 
E v > 5x 10 7 GeV, the threshold corresponding to proton energies of 10 9 GeV (a daughter 
neutrino carries ~ 5% of its parent proton's energy [39]). Below this, the WB bound is 
only an extrapolation, since the cosmic-ray flux below 10 9 GeV is totally dominated by 
the galactic component. Second, by its definition the WB bound is applicable only to 
optically thin sources. If neutrino-emitting opaque objects existed, their contributions 
might sum to produce a neutrino flux exceeding the WB bound. Potential sources 
include baryon-rich GRBs [42] and starburst galaxies [43]. 

We also adopt a model for high-energy neutrino production via shocks in GRBs [41]. 
In contrast to the WB bound, it is a prediction; see the red dashed line in Figure 1. 

2.2. Atmospheric conventional neutrinos 

While the astrophysical neutrino flux is totally unknown, there is a guaranteed and well- 
measured neutrino component — atmospheric neutrinos. These arise from the decays of 
mesons produced by cosmic rays striking the upper atmosphere. Neutrinos coming 
from pion and kaon decays form the 'conventional' component, which is well-studied 
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Figure 1. The differential flux of high-energy neutrinos (plus anti- neutrinos) 
considered in the present study, summed over neutrino flavors. The shaded gray, 
green and yellow regions indicate the uncertainty ranges for the fluxes of conventional 
atmospheric, prompt-decay atmospheric, and astrophysical neutrinos, respectively. 
The reference extra-galactic neutrino fluxes refer to the WB limit [39] (red solid line) 
and to their prediction for GRBs [41] (red dashed line). 



both theoretically [44] and experimentally [40]. Although there is no detection of 
any neutrinos for E v > 10 5 GeV, the energy range we are mainly interested in, this 
component should quite easily be extrapolated using measured data at lower energies, 
thus providing guaranteed seeds for CHAMP production. We use the model of Ref. [45] 
for the conventional atmospheric flux. As shown in Figure 1, the well-known spectrum 
of these neutrinos falls rather steeply with increasing energy. 



2.3. Atmospheric prompt-decay neutrinos 

Hadronic cosmic ray interactions with Earth's atmosphere also produce short-lived 
charmed mesons [46-54]. Even though branching ratios into these final states are 
not large, the neutrino spectrum from the subsequent decays of charmed mesons is 
quite hard as they immediately decay, before losing energy. As a consequence, the 
contribution from this 'prompt-decay' component to the total flux of the atmospheric 
neutrinos falls less rapidly with energy than the conventional component. The absolute 
normalization of this prompt flux is however still unknown, and there is a large range 
in model predictions. The main source of uncertainty is the proper treatment of next- 
to-leading order charm meson production cross sections, which strongly depend on the 
behavior of the nucleon parton distribution functions (PDFs). In our study, we consider 
a range between a smaller prompt flux from Ref. [52] and a larger prompt flux obtained 
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by using the shape presented in Ref. [53] and normalizing it to IceCube's experimental 
upper limit [54]. Such a large flux, just allowed by the data, is in fact characteristic of 
the largest model predictions among Refs. [46-54]. We show these fluxes in Figure 1. 

3. Interaction cross section 

Several TeV-scale extensions of the Standard Model feature a meta-stable massive 
charged particle. Perhaps the best-motivated scenario from a theoretical standpoint 
is supersymmetry, which provides several examples. If the LSP is very weakly 
interacting — e.g., a gravitino or a right handed sneutrino, where the interactions with 
the rest of the supersymmetric partners are suppressed by gravitational couplings or 
a gauge symmetry — the NLSP is generally meta-stable. Specifically, a charged NLSP 
can occur in supergravity theories with a gravitino LSP [16,55,56], gauge mediated 
supersymmetry breaking setups [57, 58] , scenarios featuring a stau-neutralino near- 
degeneracy (particularly in the so-called co-annihilation region [27, 59] ; here one can have 
a neutralino LSP), or supergravity scenarios with a right-handed sneutrino LSP [60]. 

Another TeV-scale new physics setup that naturally encompasses a meta-stable 
NLP is that of universal extra-dimensions (UED) with a Kaluza-Klein (KK) graviton 
as the lowest mass eigenstate in the KK tower [6]. In the minimal UED setup, the 
next-to-lightest KK state is usually neutral, and corresponds to the KK first excitation 
of the U(l) gauge boson, if the Higgs mass is below ^200 GeV, for any value 

of the compactification inverse radius R" 1 . However, as pointed out in Ref. [61], even 
in the minimal UED setup the next-to-lightest KK particle (NLKP) can be charged, 
and the LKP can be the KK graviton if Rr 1 < 809 GeV and m h > 250 GeV. In this 
case, the NLKP corresponds to the KK first charged Higgs mode. Alternatively, the 
boundary conditions at the orbifold fixed points can alter the spectrum, and give rise to 
scenarios with a KK lepton as the NLKP, and again, a KK graviton LKP [6] . In general, 
electroweak precision observables constrain the scale R^ 1 where the first excitation of 
a 5-dimensional UED scenario might be expected to a few hundred GeV, depending 
upon the value of the Higgs mass [62]. The analysis outlined below applies, with the 
appropriate production cross sections and energy losses, to such UED models and to 
any other similar framework featuring a meta-stable charged particle. 

We choose to work with two of the well-known and well-motivated supersymmetric 
frameworks mentioned above: gauge mediated supersymmetry breaking (GMSB) [57, 
58], and minimal supergravity (mSUGRA) with a gravitino LSP (see, e.g., Refs. [16,55, 
56]). For each framework we examine two models. One is a 'supersymmetric benchmark' 
(the SPS7 point of Ref. [63] for GMSB, and the e model of Ref. [64] for mSUGRA 
with gravitino LSP), and the other is a variant with a lighter spectrum (models I and 
II). These are essentially rescalings of the first two. By adopting benchmark models, 
not only do we consider phenomenologically viable and theoretically soundly-motivated 
setups, but we also make it easier to compare the detection technique discussed here 
with a wealth of existing phenomenological analyses of the same models (see, e.g., 
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Refs. [63,64]). As will become apparent below, for the present analysis the details of 
the spectrum of the heavy supersymmetric particle pair produced is crucial. Assuming 
a degenerate sfermion spectrum, while potentially useful to get an understanding of the 
role of the mass scale in the expected size of the signal, can potentially be a misleading 
over-simplification. 
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Table 1. The input parameters for the mSUGRA (upper pair) and GMSB (lower 
pair) models used in the present analysis and in Figure 2. 
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Table 2. The masses of the lightest stau, first-generation squark, and lightest 
neutralino and chargino, for the four models of Table 1. 

We specify the mSUGRA and GMSB input parameters in Table 1. Notice in 
particular that Model II lies on the SPS7 slope defined in Ref. [63]. In Table 2 we detail 
the four models' relevant particles masses. 

The gravitino mass need not be specified as long as the stau decay length, ct?, is 
larger than or of the same order as the Earth radius, i? e . This implies a lower limit on 
the gravitino mass mg [65], 




Rearranged, the formula implies, for instance, that for a 100 GeV stau the gravitino 
mass can be as light as 1 MeV, and for a 1 TeV stau mg has to be larger than 0.3 GeV. 
Notice that Equation (1) does not take into account the relativistic boost factor, 7, 
which is sizable for staus produced in very high energy neutrino-nucleon interactions. 
It is therefore a conservative constraint on the gravitino mass. 

In addition, the stau lifetime should be short enough to be consistent with 
limits obtained from the effects on light element abundances processed in big-bang 
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nucleosynthesis [66] (see also Refs. [67]) and from excessive distortions to the cosmic 
microwave background spectrum [68]. Of particular relevance are constraints resulting 
from overproduction of 6 Li and 7 Li [28, 69, 70], induced by catalytic effects produced by 
bound states consisting of light nuclei and the meta-stable CHAMP (here, the lightest 
stau). Notice, however, that the details of the estimate of the amount of primordially 
synthesized lithium are still under debate [71]. Also, the constraints depend upon the 
fraction of electromagnetic energy released in the decay. Conservatively, if one requires 
the lifetime of the charged meta-stable species to be shorter than 10 3 -10 4 s, as implied 
by the analysis of Ref . [72] , the gravitino mass is constrained to be approximately below 
1 GeV for a 100 GeV stau, and below 100 GeV for a 1 TeV stau. This evidently leaves 
a very wide window, of almost three orders of magnitude, for the viable gravitino mass 
range. 

To calculate the stau frux,|| we first need to calculate the stau production cross 
section as a function of incoming neutrino energy. We do this using the susy-madevent 
package [73,74], which calculates the differential or total cross section for any 2 — > n 
scattering process in the MSSM given an SLHA-conforming (standardized format for 
spectra) model input file [75]. We generate model input files using the SUSPECT 
spectrum generator package [76], which also automatically checks the generated model 
against various known precision data constraints, such as b — > s'y. Our four models do 
not conflict with any known constraints. The processes contributing to stau production 
mainly stem from tree-level u and t exchange diagrams with a slepton (l and v) plus a 

squark (u and d) in the final state, through neutralino (x°) or chargino (x ± ) exchange, as 

x ± ~ ~ x° ~ ~ ~ 

discussed in Ref. [33]. More specifically, they are: u+u(d) —> l+d(u), u+u(d) — > u+u(d), 

where x° an d X^ indicate neutralinos and charginos, respectively. The squarks and 

sleptons produced in the final state promptly cascade decay into meta-stable staus. 

We calculate total cross sections, summing over neutrino and anti-neutrino inelastic 
scattering on protons using exact matrix elements to supersymmetric pair final states for 
both charged current (CC) and neutral current (NC), employing CTEQ6L1 PDFs [77]. 
The possible final states are (anti-)sneutrino or (anti-)slepton, plus (anti-)squark. Our 
calculations are leading order, as there are no available NLO QCD corrections for 
neutrino-proton (up) scattering. Judging from the known results for proton-proton 
(pp) scattering, however, we probably make an underestimate of the rate on the order 
of 50%. In this regard our calculation is conservative. For our purposes at the relevant 
energies, the neutrino-neutron cross section is sufficiently close to the neutrino-proton 
cross section that we can ignore calculating it separately. 

An interesting side observation is that the up and up cross sections are not equal 
except at very large y/s, where low- a; quarks dominate the PDFs and are approximately 
egalitarian [78]; we reproduce this observation here. As \/s gets within a couple orders 
of magnitude above threshold, however, up dominates because the larger CC process 

|| Here the word 'stau' denotes staus and anti-staus collectively. This is because both could be produced 
by the same interaction, but cannot be distinguished at neutrino telescopes. 
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picks out valence quarks, where u dominates slightly over d. Closer to threshold this 
remains true for the CC process, but for NC vp dominates because of chirality-selection 
for the final-state squarks. Near threshold these two diverging components accidentally 
roughly cancel, resulting in o vp w o~ Pp once again. For the present purposes, however, 
we are concerned only with total rates, and ignore charge-separated subsamples, which 
would vary somewhat as a function of stau energy. 

The left panel of Figure 2 illustrates our results for the neutrino-proton scattering 
cross section a up (E u ) into any supersymmetric particle pairs, for the four models listed 
in Table 1, as a function of the incident neutrino energy. The general trends in a up (E u ) 
are consistent with those found in other analyses, see e.g., Ref. [33]: the steep rise 
in the low-energy end reflects the strong kinematic suppression associated with the 
squark-slepton pair production threshold. The subsequent rise of the cross section with 
the incoming neutrino energy depends upon the small-rr behavior of the PDFs. We 
give in Section 5, where we discuss the role of the specific supersymmetric particle 
spectrum in the determination of a vp {E u ) — > SUSY, an analytical interpretation of 
the specific power-law behavior that emerges from the numerical computation. Notice 
that compared to the optimistic toy model used in Ref. [33], theoretically-motivated 
(optimistic) supersymmetric setups appear to give a maximal a vp that is roughly one 
order of magnitude smaller in the asymptotic high-energy regime. 

In the right panel of Figure 2, we show the ratio of the neutrino-nucleon cross 
section into supersymmetric particle pairs over the total neutrino-nucleon cross section 
(that is, as apparent from the figure, always close to the purely Standard Model cross 
section). As we explain in Section 4.2, this is the physical quantity of interest, in the 
limit of Earth as a thick target for high energy neutrinos, for the computation of the 
stau flux. Beyond threshold effects, we point out that in the energy range of interest 
the branching ratio of neutrino-nucleon interactions into supersymmetric particle pairs 
lies between 10~ 4 and 1CT 3 . 

4. Flux of long-lived staus 

We devote this section to a detailed analytical treatment of the computation of the flux 
of staus produced by neutrino-nucleon collisions that might be detected at a neutrino 
telescope such as IceCube. We start, in Section 4.1, with a derivation of the differential 
flux of staus from first principles, leading to the result presented in Equation (6). In the 
following Section 4.2 we assume that Earth is opaque to neutrinos at energies relevant 
here. In this thick target approximation, we analytically show that the flux of staus 
at the detector can be computed as a simple integral, over incident neutrino energies, 
shown in Equation (10), of the product of three factors: 

(i) the differential flux of incident neutrinos (shown in Figure 1), 

(ii) the ratio of the neutrino-nucleon cross section into supersymmetric particles over 
the total neutrino-nucleon cross section (shown in Figure 2, right), and 
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Figure 2. Left: The neutrino-proton cross sections, as a function of the incident 
neutrino energy for the production of stau pairs in the four models under consideration 
here (see text and tables 1 and 2 for details on the specific models). Right: The ratio 
BR(jaN — + SUSY) of the neutrino-proton cross section into supersymmetric particle 
pairs over the total neutrino-proton cross section, as a function of the incident neutrino 
energy, for the same four models as in the left panel. 

(iii) a 'geometric efficiency' factor, to be defined below and explicitly shown in the right 
panel of Figure 3. 

When a neutrino interaction occurs, the branching ratio for stau production 
among the final states is given by the above cross section ratio. In the thick-target 
approximation, all incoming neutrinos below the horizon will interact in Earth. Since 
the staus are collinear with the incoming neutrino direction, and are produced with a 
sizable fraction of the neutrino energy, this leads to a simple but important result about 
the stau flux. Neglecting stau energy losses in matter for a moment, we see that Earth 
acts as a neutrino-to-stau converter, with a probability that is independent of direction 
(below the horizon). This gives an upper bound on the stau flux through IceCube. The 
relevant range of energies will be between threshold (~ 10 6 GeV; see Figure 2) and the 
point at which the stau flux becomes too small (~ 10 7 GeV; note the product of Figures 
1 and 2). Taking stau energy losses into account will only reduce the stau flux at the 
detector. Since staus cannot reach the detector from too far away, this means that only 
a limited range of nadir angles will be relevant, and this defines our geometric efficiency. 

With our approximations, we give a preliminary assessment of the total expected 
stau flux at the detector (Table 3). Finally, we compute the actual accurate stau flux 
resulting from the full-glory integration of Equation (6) in Section 4.3. We provide 
numerical estimates of the integrated flux in Table 4 (that fall within a factor 2 of the 
approximate results anticipated in Table 3), as well as the actual differential flux of staus 
from the different primary incident neutrino sources. Finally, we comment in Section 4.4 
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on the flux of staus predicted from nucleon-nucleon reactions. 



4-1. Formulation 

In the framework we consider here, all supersymmetric particles produced in neutrino- 
nucleon interactions promptly decay into the NLSPs — here, the stau, which is meta- 
stable: long-lived enough to propagate through Earth (with energy loss) and reach the 
detector. Our objective is to calculate the spectrum of staus after this energy loss, 
following similar principles for the spectra of neutrino- induced muons [79,80]. 
The stau electromagnetic energy loss rate is given by [32, 35] 
dE 

= a s + (3 S E S , (2) 

where X = p x / is the column depth of matter in units of g-cm -2 (i.e., density times 
distance), and the a s and j3 s terms represent ionization and radiation losses, respectively. 
We neglect discrete scattering by weak interactions, as the effect is small at the energies 
we focus on near threshold [36] . Hereafter, we use the subscript s to indicate quantities 
referring to staus. Our coordinate system locates the detector at X = and particles 
are produced at X > 0, so that the energy E s is a growing function of X (dE s /dX is 
positive). For the density profile of Earth, we use the model given in Ref. [81]. 

The ionization coefficient for staus, a s , is approximately the same as that of 
muons [35]; specifically, we use a s = 2 x 1CT 3 GeV cm 2 g _1 . Radiative losses, on 
the other hand, depend on particle mass, and we take the corresponding coefficient to 
be given by [3 S = 4.2 x 10~ 6 (m M /m s ) cm 2 g _1 [35]. By integrating Equation (2), the 
'distance' Xif = Xi — Xf traversed by the stau while its energy decreases from E l s to 
El reads 

' ft \a, + l3,Et) 
We assume that each stau produced in up interactions carries a large fraction of the 
parent neutrino energy. We denote this as E l s = (1 — y)E v and assume y = 0.5, 
independent of neutrino energy. (We discuss in Section 4.2 below the dependence of 
the stau flux on the parameter y.) At these high laboratory energies, the staus may 
be taken to be collinear with the original neutrino direction. The differential flux of 
produced staus per energy E l s and distance X reads 

(Am ax — X)a tot (E u 



t d*¥ s a slJSY (E v )d<S> u (E u ) 

- dEldx - (1 _ y)mp dEu ex P 



(4) 



where E v = E l s /(1 — y). The exponential factor takes into account the neutrino 
attenuation, mainly due to the Standard Model interactions (<7 to t ~ o'sm] see Figure 2, 
right); X max represents the value of X corresponding to the surface of Earth, which 
depends on the direction. The overall (1— y)^ 1 factor comes from the change of variables 
from the spectrum of neutrinos to that of produced staus, i.e., dE^/dEl = (1 — y)^ 1 . 
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The spectrum of staus at the detector is given by a double integral over all 
production positions and energies, subject to the constraint of having the stau energy 
be between E[ and E[ + dE{ at X = 0: 



where the function f(E l s ,X) is defined by the energy loss, Equation (3). We use the 
energy constraint to perform the E\ integration, and as a result we obtain the following 
expression: 



where E v and E\ inside the integral have to be evaluated according to the chosen E{ 
on the left-hand side and the X at that step inside the integral. From the energy-loss 
equation, Equation (3), we have 



where, again from the kinematic definition, E u = E\/(l — y). 

It is convenient to change the integration variable by dividing the differential dX 
by X and multiplying the integrand by X. Then the integration steps are in lnX, 
and in Figure 4 we show this new integrand for different nadir angles. In the left-hand 
panel, the stau energy losses in matter are neglected, so that the neutrino interaction 
and geometric effects are shown clearly. For each nadir angle (90° is at the horizon, and 
0° is along the diameter of Earth), the sharp edges at large X occur due to the boundary 
of Earth. The peaks arise because the neutrino interactions occur logarithmically near 
X = X max . Substantial attenuation of the neutrino flux is only seen at small nadir 
angles; in units of the axis, the exponential scale height is ~ 0.1. In the right-hand 
panel, the stau energy losses are now included. As expected, this prevents staus from 
arriving at the detector from too large of distances (beyond ~ 0.01-0.05 in units of the 
axis, depending on nadir angle due to the radial variation of the density profile). In this 
panel, the visual area under each curve shows its relative importance to the total stau 
flux through the detector. 

4-2. Thick target approximation 

For neutrino energies relevant for our purposes, E v > 10 6 GeV, Earth is opaque to 
neutrinos at the most important nadir angles. Thus, to a good approximation, we 
can assume that staus are produced from neutrino interactions logarithmically near 
Earth's surface. In this case, the exponential factor in equation (6) is sharply peaked at 
X = X max , i.e., exp[-(X max - X)a SM /m p ]/( y m p /o- S M) - S(X - X max ). For our change 
of variables, the integrand is proportional to 5(\nX — lnX max ) = X5(X — X max ). 








(7) 
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Figure 3. Left: Minimum neutrino energy, -E„ ;m i n (#), below which the produced staus 
cannot reach the detector, as a function of the nadir angle 9. Right: The geometric 
efficiency factor e goo (£"„), as a function of the incident neutrino energy (see Equation 
(10)). 



We then obtain 

'&E 1 , El V dEj a SM KdE>,\ dKJ "sm K 

where, in the second equality, we used the relation dE l s /dE{ = e /3Xmax from Equation 
(7) with X = X max . 

We require the staus to be more energetic than a given detector threshold E s ^. 
If the staus are relativistic enough to emit Cerenkov light, they will be detected. We 
adopt £7 Sj th = 300 GeV, since the typical stau mass we consider is ~100 GeV, but we 
note that our results change negligibly for higher or lower thresholds. For instance, a 
shift by one order of magnitude, E s ^ = 3 TeV, affects the final result by only 0.2%. 

To satisfy the detection requirement E[ > E s ^, the initial stau and corresponding 
neutrino energies must be larger than some minima, E l smin and E u ^ min , given by 

El, mm (9) = (E s , th + |) e^«) - |, E u , mm (9) = (9) 

as evident from Equation (7). We plot this minimal neutrino energy as a function of the 
nadir angle 9 in the left panel of Figure 3 for the various supersymmetry models under 
consideration. Note that the crucial quantity here is the stau mass, hence we obtain 
the same result for models I and II. In the thick target approximation, the flux of staus 
reaching the detector from below is thus given by 

r r 00 .df^ 

F„= dn dEi 



e<n/2 JE Sith ° dEi 
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Model 


WB bound 


WB GRB 


Atm. Prompt 


Atm. 


Conv. 


I 


< 3.2 


0.20 


0.012-0.73 


0.0023 


-0.0038 


II 


< 1.2 


0.066 


0.0028-0.16 


0.00028- 


-0.00048 


SPS7 


< 0.88 


0.045 


0.0018-0.099 


0.00014- 


-0.00024 


e 


< 0.71 


0.034 


0.0011-0.062 


0.000069 


-0.00012 



Table 3. Stau fluxes from various neutrino sources for the four benchmark 
supersymmetry models, in the thick target approximation, in units of km -2 yr _1 . 

= / dil r AE U (^) (?*™L) {E v - 
Jd<n/2 Jo V di W V °sm / 

■-"f^CtJte)^ (10) 

where 6 is the step function, and in the second equality, we have simply used definitions 
given above. In the last equality, we are defining a 'geometric efficiency 1 factor 
e geo (E u ), assuming that the incident neutrino intensity § u is isotropic. We show this 
geometric efficiency factor in the right panel of Figure 3. Thus, under the thick target 
approximation, the detection flux of staus can be divided into three independent factors: 

• incident neutrino spectrum d§ u /dE u 

• cross section ratio o"susy/csm 

• geometric efficiency e geo 

These three factors are illustrated in figures 1-3. After integrating over the neutrino 
energy in Equation (10), we obtain approximate stau fluxes at the detector, summarized 
in Table 3 for the ranges of incident neutrino fluxes given in Section 2. 

So far we worked with the assumption that y = 0.5; i.e., each stau carries half of 
the incident neutrino energy. However, especially when including the complex pattern 
of chain decay of squarks and sleptons into staus, a smaller fraction of the maximally 
available energy is expected to be carried by the staus, implying a larger value for y. 
The precise value for y is model-dependent, and its detailed evaluation is beyond the 
scope of the present analysis. We thus investigate the dependence of the stau flux on 
the y parameter, for simplicity under the thick target approximation. Our calculations 
show that the flux of staus, with an incoming flux saturating the WB bound, and with 
our benchmark supersymmetric model I, is 3.2, 2.7, and 1.5 km -2 yr -1 for y = 0.5, 0.7, 
and 0.95, respectively. This shows that the stau flux changes at most a factor of ~2 for 
a wide range of y, which is well within other model uncertainties. The weak dependence 
we find stems from the fact that a larger value of y requires larger neutrino energies 
to produce staus with a certain energy. In turn, at larger energies the incident flux is 
smaller, but the cross section ratio is larger. The cancellation of these two effects results 
in the mild dependence we find. The same argument also applies when we evaluate the 
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Figure 4. The contribution to the stau events from a given column depth X, at an 
energy E[ = 10 6 GeV, for an incident neutrino flux saturating the WB bound. Each 
curve corresponds to a different value of the nadir angle 8 as labeled. Left panel is 
shown with the assumption that there is no energy loss for staus, just for illustration 
purpose. In the right panel, proper energy losses are included. 

flux more accurately in the next subsection, where we thus again use the assumption 
y = 0.5, for simplicity. 

4-3. Results of numerical integration: stau spectrum 

We now solve Equation (6) numerically to obtain a more precise spectrum and rate 
estimate of staus at the detector, as well as to check that the approximation made in 
the previous subsection is reasonable. Before giving the final flux estimates, we start by 
investigating the generic structure of the integrand of Equation (6). 

Figure 4 shows the integrand of Equation (6) as a function of the column depth X 
for various values of the nadir angle 8, assuming E* g = 10 6 GeV. For definiteness, we 
show here the result only for a neutrino injection model saturating the WB upper limit. 
From the right panel of this figure, one can see that the stau events are dominated by 
directions with large enough nadir angle, specifically with 9 > 70°, so that the staus can 
reach the detector. On the other hand, for very large nadir angles, the contributions 
to the total event rate are modest, because Earth is not completely opaque in those 
directions, even at these high energies. 

Figure 5 shows the differential stau fluxes as a function of final stau energy. First, 
we note that nearly all of the detectable staus are well above the threshold required to 
be relativistic, which is a stau energy comparable to the stau mass. Only relativistic 
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charged particles produce the Cerenkov light that IceCube can measure. 

The energy loss of relativistic particles may be dominated by ionization or radiation, 
depending on whether the a term or the /3 term dominates in Equation (2), respectively. 
This transition for muons occurs at an energy a/f3 ~ 500 GeV. For staus, it occurs at 
an energy a factor ~ m s /m^ higher, i.e., at least 10 6 GeV. The energy loss associated 
with Cerenkov radiation is always negligible; on the other hand, the Cerenkov radiation 
per unit length is the same for all relativistic particles. Thus, for particles at any energy 
in the relativistic ionization-dominated regime, all tracks will look the same in IceCube. 

At higher energies, in the radiation dominated regime, there is additional Cerenkov 
radiation arising from relativistic electrons and positrons created in hard radiative 
processes. In this regime, one can indeed tell the energy of the primary particle by the 
intensity of the total Cerenkov radiation. For most of the relevant final stau energies 
shown in Figure 1, the staus will at most be only slightly in the radiative regime, and 
so all stau tracks going through IceCube will be indistinguishable from each other (and 
from low-energy muons). While the total energy deposited in the detector is much 
smaller for staus than it is for low-energy muons, this is irrelevant for IceCube, which 
detects only the Cerenkov light. 

We argue that low-energy but relativistic stau pairs could also be detected (albeit 
without energy measurement), giving a sizable event rate. Recall that while it might 
be difficult to distinguish between staus and muons on the basis of a single-particle 
detection, it would still be possible if we use dual-track events: since staus propagate over 
much longer distances in Earth than muons, tracks entering the detector simultaneously 
are expected to be well-separated [32-34]. The careful analysis of Ref. [34] shows that 
the separation distribution of stau pairs ranges widely from 50 m to 1 km, and that it 
peaks around ~500 m. On the other hand, the separation distribution for di-muons— 
the main background for the stau pair track — peaks at around 10 m, and essentially no 
di-muon events with > 50 m separation are expected. Therefore, this criterion rejects 
almost all background di-muon events, but would capture a large fraction of stau events 
(typically > 50%). Thus, we are interested in the stau flux integrated over energies 
larger than the relativistic threshold — which is very small, -E^th — 300 GeV — and where 
our results are least sensitive. This is clear from Figure 5. 

Table 4 shows the expected stau flux at the detector, obtained by integrating the 
spectrum above 300 GeV. Comparing with the results of Table 3, where we used the 
thick target approximation, we find the two results are consistent with each other within 
a factor of 2, justifying the reliability of the thick target approximation. 

Note that the flux of conventional atmospheric neutrinos is not totally isotropic, 
but peaks in the horizontal direction by almost one order of magnitude compared to 
other directions [45]. (The prompt flux is isotropic.) This is an important effect because 
most of the staus reaching the detector arose from neutrinos from horizontal directions 
(Figure 4). Therefore, our results in Tables 3 and 4 for atmospheric neutrinos would 
be larger by a factor of ~3, as these results were obtained with a direction- averaged 
incident neutrino flux. 
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Figure 5. The differential stau flux at the detector for the four benchmark 
supersymmetry models. Each line corresponds to a different incident neutrino flux 
model: WB upper limit (blue solid), WB GRB (red dashed), atmospheric prompt 
(green dot-dashed), and atmospheric conventional (magenta dotted). 



Model 


WB bound 


WB GRB 


Atm. Prompt 


Atm. 


Conv. 


I 


< 2.2 


0.13 


0.0070-0.41 


0.0010- 


-0.0018 


II 


< 0.93 


0.049 


0.0019-0.11 


0.00017- 


-0.00029 


SPS7 


< 0.70 


0.035 


0.0012-0.070 


0.000091 


-0.00015 


e 


< 0.57 


0.027 


0.00080-0.047 


0.000049- 


-0.000082 



Table 4. The stau flux from various neutrino sources for the four supersymmetry 
benchmark models, obtained from the exact numerical integration of Equation ( 6), in 
units of km~ 2 yr _1 . 



As a consequence, given that the predicted stau flux could be as large as ~1 km -2 
yr _1 and that one expects essentially no background from muon pair events, the search 
for stau pair tracks is warranted in the actual data. Our discussion here shows that a 
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significant fraction of stau events could possibly come from atmospheric prompt-decay 
neutrinos (if the flux is close to the current upper bound), regardless of the assumed 
supersymmetry models. 

4-4- Stau flux from nucleon-nucleon collisions 

The stau pair flux <3> s produced from a differential flux of primary high-energy nucleons^f 
d$jv/d-EAr colliding with atmospheric nuclei is given by 
d$ s <taw,susy d$Ar dE N 

— ■> (11) 

dE s o~NN,tot dE]y dE s 

because the thick target approximation (introduced in Section 4.2) is very good for 
NN interaction. The symbols a nn, tot and ctjvjv.susy indicate the total nucleon-nucleon 
interaction cross section and the cross section into any super symmetric particle pair, 
respectively. To reiterate, since direct decays into gravitinos are strongly suppressed by 
gravitational couplings, all final state R parity odd particles decay into the NLSP, i.e., 
lightest stau pairs. For the total NN cross section, we assume the parameterization for 
the hadron-air total cross section [82] 

^, 185 + 1 ,3 log (^) + 0.08V(^), (12) 

where A ~ 14.6 is the average number of nucleons in a nucleus of air. We approximate 
the nucleon-nucleon cross section into supersymmetric particles with the proton-proton 
cross section, and we compute the latter using Prospino2.0 [83]. For the incoming 
nucleon flux we use the estimate in figure 1 of Ref . [82] . 

Quark-anti- quark processes can produce stau pairs directly, unlike neutrino-quark 
processes, where the final state has a larger threshold as the final state must contain a 
typically-heavier squark. Hence the kinematic threshold, as a function of the primary 
particle energy, is lower in nucleon-nucleon collisions than in neutrino-nucleon collisions. 
Also, the subsequent occurrence of various supersymmetric particle thresholds at larger 
and larger masses, including particles featuring large degeneracy factors (such as 
squarks), implies a more rapidly growing behavior for the ojvat.susy cross section than 
that for 0Vtv 5 sijsy- As a consequence, we find that for the models under consideration 
here, E p dN s /dE p is almost constant over several orders of magnitude in E p . 

The final flux of staus is, however, dramatically suppressed by the ratio 
o 'n n, susy /&nn, tot [84], even taking into account multiplicity effects in stau pair 
production or proton re-interactions. In particular, for the models we consider here, 
where the strongly interacting supersymmetric particles are typically much more massive 
than the NLSP, the combination of threshold effects and of the rapidly decreasing flux 
of incident primary protons leads to dramatically less optimistic predictions than those 
recently reported in Ref. [37]. There, the authors considered squarks and gluinos with 
extremely low masses (150 and 300 GeV), while the theoretically-motivated benchmark 

% Secondary nucleons and other hadrons contribute a small fraction of the incoming flux, and taking 
them into account does not affect our conclusions. 
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models we use here feature squark and gluino masses between 600 GeV and 1 TeV. 
While we agree with the numerical results reported in Ref. [37] when making the 
same assumptions on squark and gluino masses, we obtain much lower figures for the 
benchmark models we adopt here. Namely, we find that for the two most optimistic 
models, I and II, we predict a stau flux in IceCube from nucleon-nucleon interactions 
of 10~ 4 and 3 x 10~ 6 per km 2 per year, respectively — much lower than even the 
contribution from conventional atmospheric neutrinos. We believe that this relative 
smallness compared with that from the atmospheric incident neutrinos would be a rather 
model-independent feature. 

5. Role of the supersymmetric particle spectrum 

In the previous sections we focused on specific supersymmetric models. We now wish to 
address the model-independent question of how the stau pair rate at neutrino telescopes 
depends upon the supersymmetric particle masses. As we already pointed out, the stau 
pair flux depends upon the uN — > SUSY cross section: the relevant masses entering the 
cross section are those of the heavy pair produced, and those of the supersymmetric 
partners of the electroweak gauge bosons which mediate the charged- and neutral- 
current interactions responsible for slept on-squark production. In addition, the number 
of produced stau pairs as a function of incoming neutrino energy crucially depends on 
the final state kinematic threshold: the larger the latter, the smaller the flux of incoming 
neutrinos that can lead to stau pair production, hence a smaller expected stau pair rate. 

To explore quantitatively the statements above, we employ a convenient 
phenomenological parameterization of the supersymmetric setup at the low-energy scale, 
and no longer rely on a specific supersymmetry breaking framework. The left panel 
of Figure 6 shows the stau production cross section variation in high-energy neutrino- 
proton collisions in the plane defined by the slepton (x-axis) and squark (y-axis) masses. 
For definiteness, we fix the relevant gaugino (respectively, bino and wino) masses to 
m! = 1 TeV and m 2 = 2 TeV. We choose fi = 1 TeV, m^ = 500 GeV, tan/3 = 10, set 
all trilinear scalar couplings to zero, and further assume that the soft supersymmetry 
breaking scalar masses of sleptons and (separately) of squarks are degenerate, CP- 
conserving and flavor-diagonal. In the figure, we show iso-level curves at fixed values 
of the production cross section for staus at an incident neutrino energy of 10 8 GeV, 
in order to avoid production threshold effects. As the figure illustrates, cross section 
variation is very mild well above threshold. Quantitatively, the effect varies within little 
more than one order of magnitude for scalar masses varied between 100 and 1000 GeV. 

In the inset, we show how the cross section scales the masses of t-channel 
supersymmetric particles (neutralinos and charginos) exchanged in squark-slepton pair 
production from neutrino-proton collisions. Neutralino and chargino masses are entirely 
determined at tree level by the gaugino soft supersymmetry breaking masses mi and 
m 2 and the higgsino mass term /z, and by tan/3. For simplicity we assume a common 
' gaugino /higgsind mass scale, m^, defined as the common value of mi = m 2 = fi. We 
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Figure 6. Left: Isovalue curves (labeled) of the neutrino-proton cross section for 
supersymmetric particle pair production at incident E v — 10 8 GeV, in the plane defined 
by the slepton and squark masses (for simplicity we assume degenerate sleptons and 
degenerate squarks; see text for details of the supersymmetric models). The inset 
illustrates the behavior of the cross section, again at E v = 10 s GeV, as a function 
of a common gaugino/higgsino mass (mi = to 2 = /i). Right: Contributions from 
conventional, prompt decay and extra-galactic high-energy neutrinos to the total stau 
flux, as a function of the supersymmetric scalars' (degenerate) masses), normalized to 
be relative. 



employ a common scalar mass mj for both sleptons and squarks of 300 GeV, and set all 
other low-scale supersymmetric parameters as in the rest of the figure (m^ = 500 GeV, 
tan/? = 10, all trilinear scalar couplings zero). As we illustrate in the inset, beyond 
the scalar mass scale (m^ < ms) where kinematic effects play a non-trivial role, cross 
section scaling goes like the gaugino/higgsino mass scale to the power —2. This can be 
analytically understood, since 

1 ix r^^^r Ax r AQ i i* ■«.*.<?)] 



irrig/s 



x -1 / 3 s dx 



. 2 • 

m^ + xsmL 



(13) 



where, in Equation (13), we made use of the fact that in the large E v regime, 
[x ■ q(x, Q 2 )} ~ x -1 / 3 [85]. One thus gets: 

1 / s \ 1/3 

^~^U4J ' (14) 

which explains both the scaling in the inset of Figure 6 and of a up in Figure 2. 

As the supersymmetric particle spectrum gets heavier, not only does the neutrino- 
proton cross section become smaller (left panel of Figure 6), but more importantly 
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the shift in incident neutrino energy threshold for stau production strongly suppresses 
the final stau flux. This depends on the dramatic dependence of the flux of incident 
neutrinos on energy, as illustrated in our Figure 1. The right panel of Figure 6 quantifies 
this trend. There we show the relative contribution to the total stau flux from various 
incident neutrino fluxes as a function of the common squark and slepton masses ms- 
The fluxes are normalized to that resulting from the incident WB upper limit of extra- 
galactic neutrino flux and m s = 100 GeV. We set all the supersymmetric parameters 
to the same values as in the left panel. Comparing the relative flux for various origins 
summarized in Table 4 with that shown in the right panel of Figure 6, we can roughly 
estimate that the four benchmark models correspond to ms ~ 500 GeV in the context 
of this phenomenological approach. In addition, recalling that the benchmark models 
predict stau flux of ~1 km -2 yr" 1 for WB neutrino bound, this suggests that one can 
expect a stau flux well above 1 km~ 2 yr _1 , provided slepton and squark masses are 
smaller than 500 GeV. 

The different neutrino flux scaling with energy dictates that the relative importance 
of the various neutrino sources depends on the mass scale of the particles produced 
in the neutrino-proton collision. In particular, while with a very light spectrum the 
contribution from conventional atmospheric neutrinos can be comparable to (or even 
larger than) the extra-galactic neutrino component, for ms > 0.5 TeV the contribution 
from atmospheric neutrinos becomes negligible. Note that in that mass range the 
supersymmetric particles are so heavy that the overall stau flux is extremely suppressed, 
and likely undetectable. On the other hand, the figure illustrates that, in principle, 
prompt decay neutrinos can be the dominant source of staus for almost any value of 
m s if the extra-galactic neutrino flux is close to the GRB-derived range (dashed red 
line in the figure) rather than the WB upper limit. In addition, even conventional 
neutrinos contribute a stau flux of the same order of magnitude as that expected from 
the WB upper limit on astrophysical neutrinos, as long as the supersymmetric scalars 
mass scale is below 200 GeV. If both the prompt-decay neutrino flux and the extra- 
galactic neutrino flux are maximal, prompt neutrinos contribute at the same level as 
extra-galactic neutrinos for < 400 GeV, and dominate for ms below 200 GeV. Since 
a detectable signal is expected only for a light supersymmetric spectrum, this leads us 
to the following prediction: if the signal discussed here is indeed detected, a very sizable 
fraction of it will originate from conventional and prompt-decay neutrinos. 

6. Conclusions 

We reassessed the flux of meta-stable staus produced by neutrino-nucleon and nucleon- 
nucleon interactions that might be detectable at km 3 neutrino telescopes. We derived the 
flux of staus from first principles, and showed that, under the approximation that Earth 
is opaque to very high energy neutrinos, the number of staus at the detector is given by 
a simple integral over the neutrino energy of the product of three factors: the incident 
neutrino flux, the ratio of the neutrino-nucleon cross section into supersymmetric 
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particle pairs over the total neutrino-nucleon cross section, and a geometric efficiency 
factor. We showed that this approximation reproduces an exact numerical computation 
within a factor 2, which in turn is much better than the level of our knowledge of the 
first two factors entering the stau flux computation — namely, the incident neutrino flux 
and the features of the supersymmetric particle setup. 

We focused on each of the factors relevant to the computation of the final flux 
of staus. We concentrated on four well-motivated supersymmetric benchmark models, 
and independently evaluated the relevant production cross sections. We pointed out 
that previously-neglected atmospheric neutrinos from prompt charmed meson decays 
could give a potentially large stau flux, even in the absence of a (yet to be discovered) 
astrophysical high-energy neutrino flux. This will depend on the prompt atmospheric 
neutrino flux being at the upper end of the theoretically expected range. Nucleon- 
nucleon processes, even for the most optimistic benchmark models, would not contribute 
sizably to the final stau flux. Finally, we numerically and analytically studied how 
the relevant cross sections depend on the supersymmetric model mass spectrum, and 
how the relative importance of primary neutrino sources depends on the mass scale 
of supersymmetric scalars. In particular, we predict that if the signal discussed here 
is indeed detected, a very sizable fraction of it would originate from conventional and 
prompt-decay neutrinos. 
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